Charge correlations in the weakly doped t — J model calculated by projection 

technique 
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We study frequency- and wave-vector dependent charge correlations in weakly doped antiferro- 
magnets using Mori-Zwanzig projection technique. The system is described by the two-dimensional 
t-J model. The ground state is expressed within a cumulant formalism which has been successfully 
applied to study magnetic properties of the weakly doped system. Within this approach the ground 
state contains independent spin-bag quasiparticles (magnetic polarons). We present results for the 
charge-density response function and for the optical conductivity at zero temperature for different 
values of t/J. They agree well with numerical results calculated by exact diagonalization techniques. 
The density response function for intermediate and large momenta shows a broad continuum on en- 
ergy scales of order of several t whereas the optical conductivity for lo > is dominated by low 
energy excitations (at 1.5-2J). We show that these weak-doping properties can be well understood 
by transitions between excited states of spin-bag quasiparticles. 



I. INTRODUCTION 

The discovery of the high-temperature superconductiv- 
ity has increased the interest in strongly correlated elec- 
tronic systems. Many properties of the superconducting 
cuprates are still not completely understood, an effec- 
tive theory capable of modeling the ground state and 
the excitations of these materials remains an outstand- 
ing problem. Among the interesting features in the op- 
tical response of the cuprate materials are mid-infrar 
(MIR) structures at energies at 0.2-0.5 eV, see e.g. I 
Their origin has been discussed extensively and is not 
yet completely clarified. 

The one-band Hubbard model and the t-J model are 
two candidates for a description of the copper-oxide 
planes. The purpose of this paper is an analytical study 
of the dynamical charge response and the optical conduc- 
tivity in the weakly doped t-J model at finite energies and 
zero temperature. 

Theoretical progress in this field is mostly based on 
numerical techniques, especially on exact .-jdiagonaliza- 
tion methodstjB and Lanczos calculations!! However, 
the small system sizes presently accessible to numerical 
methods leave many problems unresolved. So it is often 
not obvious whether the structures observed in the spec- 
tra obtained from small-cluster calculations are finite-size 
effects or bulk properties. Finite-size scaling is hard to 
perform for two-dimensional systems because of the lim- 
ited number of cluster sizes. Small but finite momenta 
and hole concentrations cannot be investigated by nu- 
merics. 

The spin response of doped antiferromagnets has also 
been studied in a number of analytical papers but only 
few authors have investigated the optical -conductivity 
and the charge density response function3Tl3. The low- 
energy part of the optical conductivity has been found 



to fall off slower than predicted by Drude theory; a pro- 
nounced MIR structure is visible especially at zero or 
low temperature (T <C t, J) and small doping (S < 10%). 
Some of the analytical results indicajfcfs-sharp peaks in the 
density response,at large momentall'EI. A recent slave- 
boson approacnid yields one broad structure at high en- 
ergy in this regime which is consistent with numerical 
studies. Besides the above mentioned fermionic models 
also boson models have been studied, e.g. bosons (rep- 
resenting holes) coupled to a fluctuating gauge fieldEj. 
This coupling leads to an incoherent density fluctuation 
spectrum at finite temperatures. 

The-SMatem we are interested in here is the 2D t-J 
modeMH 



H = 



.) + J 



(ij) 



(8,8,-25*). 



(1) 



Sj is the local electronic spin operator and rii the electron 
number operator at site i. The symbol (ij) refers to a 
summation over pairs of nearest neighbors. At half fill- 
ing the t-J Hamiltonian reduces to the antifcrromagnetic 
Heisenberg model. The electronic creation operators cJ CT 
exclude double occupancies: 



i(i 



(2) 



The t-J model (§ is by now believed to describe the rel- 
evant low-energy degrees of freedom of the copper-oxide 
planes in the high-T c materials. So we expect our results 
to coincide with experimental observations in the under- 
doped compounds at low energies only (1 eV and below). 
Excitations at higher energies of course are not covered 
by the t-J Hamiltonian. 

We study the time- and wave- vector dependent charge- 
charge correlation function defined by 
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G pp (k,t) = (Mpi f $H)\fo) 



(3) 



The corresponding Laplace transform can be written as 

(4) 

z = lu + irj, rj — ► . 



G pp (k,c) = (V-o|Pk t ^ T ^IV'o} 

Z Lj 



gikR.i t j g 



Here, pg = E q , CT c k+ q , CT c qa = E 
Fourier-transformed charge density operator. |-0o) de- 
notes the exact ground state of the system, and z is the 
complex frequency variable. The Liouville operator L is a 
superoperator defined by LA = [H, A] _ for any operator 
A. At zero temperature the optical conductivity cr(u>) 
for frequencies u> > is related to the charge response 
function G pp (\c,u>) as follows: 



Re <j{(jj) 



Um u;lmG pp (k,Lu) 
cN k^O 



k 2 



(5) 



where N is the total particle number. Eq. (||) follows 
from the continuity relation Lp^ — k • jk where jk is the 
Fourier-transformed charge current density operator. 

Due to strong correlations usual diagrammatic tech- 
niques based on Wick's theorem can not be easily applied 
to solve the Hamiltonian (Q). Neither the first nor the 
second part is bilinear in fermion operators, and the cre- 
ation and annihiliation operators c^ CT and Ci a do not obey 
simple anticommutation relations. For this reason, non- 
standard analytical methods like variational wavefunc- 
tions, coupled-cluster methods, slave-boson and slave- 
fermion techniques or 1 /N expansions have been applied 
to the 2D t-J model and related strongly correlated sys- 
tems. ■— - ■— - 

In the following, a projection techniqueEjlill approach 
based on the introduction of cumulants to evaluate 
charge density response functions is presented. We fo- 
cus here on the finite-energy contributions: Since we do 
not consider the diamagnetic part of the current we do 
not obtain results for the response at zero energy, so the 
present calculation does not yield a Drude-like contribu- 
tion to cr(u>). In a recent letterEa we have published first 
results for G pp (k, to) and a(w) using the present method. 
Here we employ an improved set of projection variables 
and present more details of the calculation. At inter- 
mediate and large momenta we find for the density re- 
sponse function a broad continuum of excitations on en- 
ergy scales of several t. The optical conductivity for finite 
u> is dominated by a small number of peaks at low en- 
ergies of order J. The features and the different scaling 
behavior of these spectra can be explained in terms of in- 
ternal degrees of freedom of the spin-bag quasiparticles. 

The paper is organized as follows: In Sec. II pKJariefly 
sketch the cumulant method proposed in refs.E3~LiJ and 
show how to calculate dynamical correlation functions. 
The description of the ground state of the weakly doped 
t-J model within the cumulant formalism is subject of 
Sec. III. The employed ground state wavefunction con- 
sists of independent hole quasiparticles (magnetic po- 
larons) moving on an antiferromagnetic background. The 



choice of appropriate dynamical variables for the projec- 
tion technique is shown in Sec. IV. The variables are con- 
structed from path operators which form the hole quasi- 
particles. In Sec. V we present results for the charge- 
charge correlation function and the optical conductivity. 
They will be compared with results from numerical in- 
vestigations found in the literature. A discussion of the 
results will close the paper. 



II. CUMULANT METHOD FOR CORRELATION 
FUNCTIONS 

A recently introduced approach for calculating exoeca. 
tation values and dynamical correlation functionsEjEil 
is based on the introduction of cumulants. Provided 
that the Hamiltonian of the system can be split into 
an unperturbed part Hq and a perturbation Hi, H = 
Ha + Hi, with eigenstates and eigenvalues of Hq known, 
this method uses the decomposition 



-XH 



-X(Hi+L ) 



-XH 



(6) 



This relation can be proved by comparing the equations 
of motion of cither side with respect to A. Lq is the 
Liouville operator corresponding to Hq, defined by the 
relation LqA = [ifo,-A]_ for any operator A. Let us de- 
note the ground state of the unperturbed Hamiltonian 
H by \(f> ) and its energy by e 



H \(f>o) = e \<j>o). 



(7) 



Here we are interested in calculating dynamical corre- 
lation functions of operators B v 



G vp {lo) = (^ \SB+^—SB^ |Vo) 



(8) 



where we have introduced 8B U = B v — ($>q\B v |V>o) ■ Using 
(||) one can showtll that these correlation functions can 
be rewritten as 



n+ hi 



(9) 



The operator f2 has similarity to the so-called wave oper- 
ator (or Moeller operator known from scattering theory) . 
Within cumulants it transforms the ground state |<^ ) of 
the unperturbed system into the exact ground state |^o) 
of H. Explicitly it is given by 



Q = 1 + lim 



1 



x-o x - (L + Hi) 



-Hi 



(10) 



The brackets (<^>o I I <?^o ) c denote cumulant expectation 

values formed with \4>o). The dot ■ indicates that the 
quantity inside (...)' has to be treated as a single entity 
in the cumulant formatipH. For a detailed discussion of 
cumulants see e.g. KuboE3. 
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The relation (Q) can be applied to either weakly or 
strongly correlated systems because its use is indepen- 
dent of the operator statistics , i.e., it is valid for fermions, 
bosons or spins. Note that cumulants ensure size consis- 
tency for any subsequent approximations for the wave 
operator f2. nri 

Using Mori-Zwanzig projection techniquet3il3 one can 
derive the following set of equations of motion for the 
dynamical correlation functions G vpi {ijj): 

^ (z<V - ay - £VM) G vi j,(u>) = Xw ( n ) 

V 

X-qv, ^r/f, and Y, vv are the static correlation functions, fre- 
quency terms, and self-energies, respectively. They are 
given by the following cumulant expressions: 



a.',,,. = > > | n* bI(lb x y n |0 O ) C Xxl , (12) 



A 



JV(w) = E^I Qt ^ ( LQ Z -QLQ LB >) ° l0o) 



X v ^ is the inverse matrix of Xv^i , and Q is given by 

Q = i-p, p = J2 b ^\M c {fapBt • ( 13 ) 



P denotes a projection operator projecting onto the sub- 
space of the Liouville space spanned by the operators B„, 
Q projects onto the complementary subspace. 

The described cumulant version of projection tech- 
nique has a conceptual advantage compared to stan- 
dard projection technique. Using projection technique 
the Laplace transform of a correlation function may be 
written as a continued fraction expansion (or a set of 
equations of motion) describing the dynamics of the sys- 
tem. The continued fraction contains expectation val- 
ues which are static quantities. These expectation values 
have to be evaluated with the ground state of the inter- 
acting system. In standard projection technique this has 
to be done separately using mean-field methods or per- 
turbation theory. In contrast, in the cumulant approach 
static and dynamical aspects of the system are treated 
along the same lines. For a further discussion seeEl. 



III. GROUND STATE WAVEFUNCTION 



In the following 5 denotes the hole concentration away 
from half filling. Our system with TV lattice sites pos- 
sesses M = SN dopant holes. The Hamiltonian is de- 
composed into Hq and Hi as follows: 



Hq — Rising — J J]] {SIS* 



riirij 



) + J(N - 2M) 



H 1 = H f + H A 



(14) 



J 



<ij> 



The ground state |</>o) of the unperturbed Hamiltonian 
Hq is an antiferromagnetically ordered Neel state with 
M holes. The holes have fixed momenta k m and are 
located on the sublattice a m (a m =f, |) 




(15) 



\(t>Neel) 



In Ho we have shifted the energy level so that 
(ct>o\Ho\fa) = 0. 

Within the cumulant method we employ an exponen- 
tial ansatzcil for the wave operator fl, fl = e s . Here, 
the operator S contains the effect of the perturbation 
Hi onto the unperturbed ground state \4>o)- We have to 
consider two parts of H±, the spin-flip term H± and the 
hopping term H t . 

The hole motion_pmcesses induced by H t are described 
by path operatorsptta which leads to the picture of spin- 
bag quasiparticlescJ. Here we define path concatenation 
operators A n ^ (i) where i denotes a lattice site, n the path 
length, and £ the individual path shape: A n ^(i) operat- 
ing on the Neel state with one hole at site i, c^i'Neei), 
moves the hole n steps away and creates a path or string 
£ of n spin defects attached to the transferred hole. For 
n = 1 there are 4 different path shapes (Fig. 1), for n = 2 
there are 12 and so on. Explicitly, the operators A n ^(i) 
for sites i on the f-sublattice are defined by 



A 2 , s (i) = J2*T S Ml R U' (*,i€t,j,mei) (16) 

3 1 
jtm 



Now we turn to the t-J model at weak doping. The 
description of the ground state within the cumulant for- 
malism used here has been shown to produce reasonable 
results for the doping dependence of the staggered mag- 
netization and the spin-wave spectmip of the antiferro- 
magnetic phase of the doped systcmEj. 



The operators A n> £ for sites i on the 'down' sublattice are 
defined analogously with all spins reversed. The matrices 
R-i ...it allow the hole to jump along a path of shape £: 



1 ii, ...,i n connected by pathshape£ , . 
•'> 10 otherwise ' ^ ' 
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FIG. 1. Path shapes of length 1 created by the operators Ai,£ acting on a hole at site i in the f-sublattice. 



Having defined the excitation operators A n ^(i) the 
wave operator used for the cumulant formalism takes 
the form 



n max myi 



n = cxp ^2 X] -w-'W 

\n=l {=1 



(18) 



Note the additional 
riji-)(l — ny) which 



with parameters A n ,£ yet unknown 

definiton A (i) = \ J2 a hA a = (1 
is the only "path" with zero length (£ = 1). This oper- 
ator is a projection operator on the empty state at site 
i. The path operators A n £ commute with each other 
because they only contain spin-flip operators destroying 
Neel order. Applying the operator Vl (|lS|) to the unper- 
turbed ground state ( |l5| ) adds to each hole a cloud_nLspin 
defects leading to a spin bag or magnetic polarorO'Lj. In 
the ansatz ([18]) we use separate coefficients X n ^ for all in- 
dividual paths 4tti£- This is an extension compared to 
the ansatz of ref.c3 where all paths with the same length 
n had been given the same weight. The additional de- 
grees of freedom improve the ground-state wavefunction 
for non-zero hole momentum taking into account that the 
quasiparticle state is not exactly s-like, see also Sec. V 
C. 

The transverse part H± of the magnetic exchange cre- 
ates or destroys pairs of spin defects and therefore can 
change the length of a given path of spin defects by 2. 
This process leads to a coherent motion of the hole quasi- 
particle through the latticeoO. 

Additional spin fluctuations in the antiferromagnetic 
background will be neglected here. The results show 
that this approximation gives already reliable results for 
the charge density response of the system because the 
charge dynamics is mainly determined by the motion of 
the hole quasiparticles. However, for a description of 
magnetic properties or the spin response of the system 
these gxflund-state spin fluctuations would be essential, 
see e.g£3. 

The coefficients A n> £ of the ansatz ( |l8| ) are determined 
by the condition for f2|0o) to be an eigenstate of H. Fol- 
lowing ref.E-j one arrives at a non-linear set of equations 
for the ground-state energy E and the coefficients A„^: 



E Q = <0 o |™|0o) c , 



Al^Hni^r (19) 



with f2 given by (|l^). The cumulant expectation val- 
ues have to be taken with respect to the unperturbed 
ground state (pi]). The set of equations ( p~9[ ) together 
with the additional approximation of independent hole 
quasiparticles leacLs-jto a generalized eigenvalue problem. 
For details see refo. The total ground-state energy de- 
pends on the initial momenta of the hole quasiparticles 
k m in (|l^). The energy minimum for the one-hole prob- 
lem is located at k = (±7r/2, ±7t/2) which leads to a hole- 
pocket Fermi surface. Further properties of the one-hole 
spectrum are discussed in Sec. V C. The used picture of 
independent quasiparticles is appropriate for small hole 
concentrations, i.e., for a dilute gas of holes moving in an 
antiferromagnetic background. 



IV. DYNAMICAL VARIABLES FOR CHARGE 
CORRELATIONS 



To calculate dynamical charge correlation functions we 
employ the cumulant version of Mori-Zwanzig projection 
technique as described in Sec. II. At first we have to 
choose a set of relevant operators B v . Here we are going 
to neglect the self-energy terms E„ M from ([l2|). This can 
be done using a sufficiently large set of operators B u to 
cover the charge dynamics we are interested in. 

One variable to be included is the charge den- 
sity operator itself. Note that the particle density 
= J2icr elkRi cL^i<r an d tue h°l e density p£° le = 
\ Hia elkRi Cia4a = X« CT e lkR >A {i) are equivalent quan- 
tities according to 



p^ e +Pk' = E eikR4 (^A+4^) 

= NS ki0 . (20) 

Note that J^a (|<W4<r + 4<r<W) i s tne projector on the 
empty and singly occupied states at site i. 

As first variable we therefore choose p^° le - Additional 
dynamical variables B v can be deduced from the action 
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the Liouville operator to the first variable p^° le , i.e. from 
Lp£° le , L(Lp£° le ) etc. All variables should provide par- 
ticle number conservation because H and the density op- 
erator do not change the total number of particles. 

The hopping part of L creates strings of spin defects 
when applied to p^° le , i.e., L™ p^° le is a sum of path con- 
catenation operators up to length n. The spin- flip part 
L± creates or destroy pairs of spin defects. It changes the 
length of a spin-defect path by 2. Additional processes 
connecting paths of different holes, i.e. hole-hole inter- 
actions, will be neglected consistently with the ground 
state evaluation procedure. Assuming that the charge 
dynamics is mainly carried by the spin-bag quasiparticles 
we define a set of variables which contains the processes 
described above: 

i,a 

= E eikR ' B «^W (21) 

where A n ^(i) is a path operator of length n with the path 
shape £ acting on a hole at site i, as defined in the previ- 
ous section. The indices n, £, A replace the general index 
v of the dynamical variables B v in Eqs. ((8||). The opera- 
tor -B n ,£,A couples to a hole at site i destroying it, creates 
a hole at site Ri + Ra and adds the path A n £. Here, Ra 
is a vector connecting two sites on the same sublattice, 
i + A is a short-hand notation for R; + Ra- As an exam- 
ple, Fig. 2 shows the effect of the variable £?2,4.o applied 
to a path state A^^Ci^tjiNeei) ■ Variables with A/0 arise 
from the application of H± at sites at the beginning of a 
spin-defect path. Such a process shortens the path by 2 
and moves the initial site of the quasiparticle (the path 
starting point) two sites away (|Ra| = 2). The operators 
-Bn.f.A provide a coupling between the ground state and 
excited states of the spin-bag quasiparticles, see Sec. V 
C. So the relevant part of the Liouville space spanned by 
the operators -B n ,£,A contains all individual path states 
of the holes. Note that the first of these operators is the 
hole density operator itself: 

fl ,i,o(k) = = ij> ikR * c ia cl (22) 

The quantity we are interested in is therefore the diago- 
nal correlation function G vv (\s., z) with v = (0, 1, 0). 

Additional projection variables could be found by ap- 
plying the spin-flip term H± to -B n .£,A at sites along 
the path. Those variables would provide a coupling to 
states with spin fluctuations separated from the spin-bag 
quasiparticles. Such processes would describe scattering 
of the quasiparticle states with spin waves, they will be 
neglected here. 

The above set of dynamical variables-is an extension 
of the one used in previous calculationstS There we had 
included only variables with A = 0. The new variables 
with A 7^ provide additional degrees of freedom be- 
ing useful if more than one hole is present in the system. 



This can be easily understood from a comparison with 
the density response of a system of free fermions: The 
correct excitations are obtained there if one splits the 
density p k = £ q 4 +q Cq = £\ e lkR *c|c 4 into a set of dy- 
namical variables {cj^_ q Cq} where q are the momenta of 
the particles present in the ground state. An equivalent 
set of variables can be obtained from linear combinations 
of the {c£ + c q } leading to e lkRi cl +A c i } which cor- 
respond to the -B n ,£,A of the present calculation. 

I t I t I 4 t I t I 

t tit ttl t 

Hill I I t I 

t t t I t t t i t 

I J I t I I I t I 

t I t J t t I t I t 

Am c iT |(p NW > B 2 . 4 ,A=o A u c iT |(p Ni;e i> 

FIG. 2. Effect of -82,4,0 on a path of length 3 created by 
A 3 ,i. 

Having chosen the set of projection variables ( pl| ) we 
can consider the equations of motion ( [LT| ) for the corre- 
lation functions of the operators -B n ,£,A- The relevant 
terms for the dynamics are the static correlation func- 
tions and the frequency terms given in ( |l2[ ) . The arising 
cumulant expectation values can be transformed to nor- 
mal expectation values according to appendix A. Calcu- 
lating the matrix elements we have neglected hole-hole 
interaction processes and geometry effects such as spiral 
paths. Within these approximations the path operators 
form an orthogonal basis set, 

(M^AnJ'M = MS nm 5 6ri , (23) 

where M is the number of holes present in the system. 
The static matrix \ from ( |l2| ) therefore becomes propor- 
tional to a unity matrix: 

^a© • (24) 

The elements of the frequency matrix from ( |l2| ) contain 
the hole motion processes. To be short, here we only 
state the results for one hopping process: 

(M^ Bl iA (L t B miVi& yn\4> y 

= t M (5„£ tm ri+l + 8 n £+l,mr))§A.& (25) 

where 5 n £. mv +i is 1 if the path A n £ is obtained from A mr) 
by one further hopping process, otherwise it is zero. For 
more details concerwing the matrix elements we refer to 
a recent publication^. 



5 



V. RESULTS 

In the numerical calculations we have included paths 
up to length 5 which gives a set of 485 path variables 
A n £. The system has a few holes (6 <C 1) in each of the 
four hole pockets (at momenta (±7r/2, ±tt/2)). It can be 
shown that in this case four different vectors A in the 
set of dynamical variables are sufficient to describe the 
charge dynamics at low doping. Note that this approx- 
imation does not take into account changes of the hole 
Fermi surface with doping, so our treatment is strictly 
valid only in the limit of S — > 0. These Fermi surface 
changes would be important for a good description of 
the low-energy intraband excitations of the hole quasi- 
particles at non-zero doping (see below). 

We have tested maximum path lengths of 2, 3, 4, and 
5. Beyond a path length of 2 we obtained no essential 
differences in the spectra except of richer structures with 
an increasing number of projection variables. Thus we 
expect only minor changes when including longer paths, 
see Sec. V C. 



A. Charge correlation function 

Results for the charge response function G pp (k,w) at 
t/J = 1, 2.5, and 5 are shown in Fig. 3. Note that we 
obtain discrete spectra because we have neglected all self- 
energies. In the figures we have introduced Lorentzians 
with a small artificial lincwidth. The main doping de- 
pendence of the spectra is found in the intensity which 



is proportional to the hole concentration. This follows 
from the assumption of independent hole motion and 
a doping-independent spin background (rigid-band pic- 
ture) . Within these approximations all matrix-elements 
in ( J12| ) are essentially proportional to the hole concentra- 
tion 6. From ([ll]) it is seen that therefore all correlation 
functions G vll are proportional to 5. 

For smaller momentum the spectral weight is mainly 
concentrated in a peak near w = 0. With increasing mo- 
mentum we find a transfer of spectral weight to higher 
energies. The figures for different values of tj J show that 
the broad structures observed in the spectra for large mo- 
menta scale with t, e.g., the maximum spectral weight in 
the charge-response function at (tt, tt) remains at energies 
of about 6t independent of t/J. This scaling behavior will 
be discussed in the following sections. 

Next we compare our data with exact diagonalization 
(ED) results. Due to the small cluster sizes no numerical 
spectra are available for small hole concentrations. Fig. 
4. shows a comparison of our data for t/J — 2.5 with 
ED data for 5 — 25% taken from ref.H. It is remarkable 
that even for this large hole concentration we observe 
a qualitatively good agreement with our spectra: Both 
results show the weight transfer to higher energies with 
increasing momentum and a broad continuum of excita- 
tions at large momenta. Differences may be either due to 
the small number of sites in the numerical calculations 
(which leads to an energy gap in all spectra and other 
finite-size effects) or due to the neglection of relaxation 
processes (which would produce finite linewidths) and 
Fermi-surface effects (which would affect the low-energy 
spectrum for non-zero doping) in the present calculation. 




G 




FIG. 4. Left panel: charge response function from Fig. 3 for t/J — 2.5. Right panel: Exact diagonalization data from refB 
for comparison. t/J = 2.5. The numerical data have been calculated for a 4 x 4 periodic cluster and a hole concentration of 
8 = 25%. 



B. Optical conductivity 

Using Eq. (g) one can deduce the optical conductivity 
from the calculated charge response function. Results for 
Recr(w) are displayed in Fig. 5. We find a few peaks at 
low energy, the positions of these main peaks scale with 
J. These features at 1.5 - 2 J were also found in numer- 
ical studies of the t-J and Hubbard models. They are 
supposed to coincide with the MIR structures observed 
in optical spectra of high-T c superconductors being lo- 
cated at 0.2-0.5 eV (for t w 0.5 eV, J « 0.15-0.2 eV), 
see e.g.B. 



Within the present calculation we do not obtain a 
Drude-like contribution to <j(uj) since eq. (||) is valid 
for non-zero frequency only, i.e., it does not include the 
diamagnetic part of the current. (Due to the neglection 
of the self-energies which excludes scattering processes of 
the quasiparticle states and the assumption of indepen- 
dent hole motion the Drude contribution would have the 
form DS(u>).) So the present calculation cannot account 
for analytical features as, e.g., the drop-off of the Drude 
peak of finite ui. 

Comparing our results with numerical data (at higher 
doping) we again observe reasonable agreement e.g. with 
the ED results from refB for S = 12.5%. 
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FIG. 5. Optical conductivity for small hole concentration and different t/J. The spectra have been calculated from the 
charge-charge correlation function using (0). 
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C. Relation to the one-particle spectrum 

For the interpretation of the calculated spectra one 
can consider the one-hole spectrum of the Hamiltonian 
within our approximation. From the diagonalization of 
the one-hole problem in the subspace of path operators 
A n £ one obtains several bands for the spin-bag quasipar- 
ticles. As is well known, the lowest band has a minimum 
at (±7r/2, ±7r/2) and a bandwidth of about 2 J. It cor- 
responds to a quasiparticle with s-like symmetry. The 
spin-bag states in the higher bands have nodes in the 
coefficients A„^ (compare ( ^8| ) ). Note, however, that 
a classification according to angular momentum (p, d, f 
and so on) is not appropriate at least for longer paths be- 
cause the path states A n ^ do not obey simple rotational 
symmetry. The total number of one-particle states is of 
course equal to the number of considered path variables. 
Increasing the maximum path length n max adds essen- 
tially many non-s-like states to the spectrum which have 
energies between the s-like states. So our results do not 
depend strongly on the maximum path length beyond 

^"max ^ 2. 

The low-energy peak at small, but finite momenta 
is caused by excitations within the lowest quasiparti- 
cle band. These excitations are treated correctly here 
only for the limit 5 — > since we have neglected changes 
of the Fermi surface with doping. The next structures 
at low energies arise from transitions between the first 
and second/third quasiparticle band, i.e., from the s-like 
groundstate to p-like states. These structures are espe- 
cially visible at small momentum transfer, i.e., in the 
optical conductivity, see the discussion below. The high- 
energy part of the spectra corresponds to excitations to 
higher bands. So our calculation supports the discussion 
given in ref.El where only the two low-lying states have 
been considered in a simplified analytical calculation for 
the main peak found in the optical spectrum. We want 
to emphasize that taking into account all string states 
(up to a truncation length) as done in the present work 
reproduces not only this main peak but the incoherent 
continuum at energies up to 8t also found in numerical 
work. So we have mapped the charge excitations for the 
t-J model (|l|) at small hole concentration to transitions 
between bands of non-interacting spin-bag quasiparticles. 
(Not taken into account are changes of the Fermi surface 
of these quasiparticles with doping and scattering of the 
quasiparticle states.) 

The different scaling behavior of G pp (k,w) (for large 
k) and <j(uj) with t/J arises from the internal structure 
of the spin-bag quasiparticles. The width of the lowest 
band as well as low excitation energies scale primarily 
with J, i.e., the excitation energies can be written as 
at + bJ with b/a ^> t/J w 2. ..5. Therefore the exci- 
tation from the ground state to the first excited state 
has an energy of order J as seen in the optical conduc- 
tivity. At large k the energies of the structures in the 
density response scale with t. This can be understood 



from the following qualitative arguments: At large k the 
operator p^° le couples the ground state to higher excited 
states of the quasiparticle. This coupling depends on t/J 
since the quasiparticle size changes with t/J. (MappiM. 
the localized hole problem onto its continuum versioncil 
one finds that the path coefficients decrease exponentially 
with a length scale proportional to (i/J) 1 / 3 .) Thus the 
spatial distance of the real-space nodes in the coefficients 
A n> { increases with t/J (for fixed n, £), in other words, 
for a fixed spatial node distance the excitation number 
n increases with t/J. Therefore with increasing t/J the 
operator p%° le (for fixed k) couples to higher quasiparti- 
cle states. This results in structures scaling with t in the 
density response function for large k. 

It is worth noting that the non-s-like states of the 
quasiparticle play an important role for the charge dy- 
namics considered here. Neglecting them, i.e. consider- 
ing all paths with the same length as ONE dynamical 
variable, would be a bad approximation. However, for 
the one-particle dynamics these states do nearly not con- 
tribute to the spectrum: If one calculates the one-hole 
spectral function in the subspace of path operators then 
only the s-like states carry considerable spectral weight. 
The reason for this behavior lies in the symmetry of the 
Hamiltonian. 



VI. CONCLUSION 

In the present paper we have studied the charge dy- 
namics in weakly doped antiferromagnets described by 
the t-J model at zero temperature. Our ansatz for the 
ground-state wave-functionO includes mobile hole quasi- 
particles. Background spin fluctuations have been ne- 
glected. Thus the ground state has antifcrromagnetic 
long-range order independent of hole doping. We have 
used this ground state together with a cumulant ver- 
sion of Mori-Zwanzig projection technique to calculate 
dynamical charge correlation functions. 

We find structures at lower energies scaling with J 
for zero and small momentum and a transfer of spec- 
tral weight to higher energies with increasing momentum. 
At intermediate and large momenta the density response 
spectrum consists of a broad continuum at energies of 
order t. The calculations do not cover a Fermi-surface 
related low-energy response (including a Drude-like peak 
in the optical conductivity). Our results for the den- 
sity response function as well as the optical conductiv- 
ity are in reasonable agreement with recent numerical 
data obtained by exact diagonalizationsfflO (see Figs. 
3-5). However, these numerical calculations have been 
done for large hole concentrations (e.g., 25%) where only 
short-range magnetic order is present. In contrast, our 
calculations are based on magnetic long-range order and 
neglect background spin fluctuations. Thus we conclude 
that the fact whether the quasiparticles move in a long- 
range ordered background or not does not have an es- 
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sential influence on the charge response of the system at 
higher energies. The dynamics is determined by the local 
antiferromagnetic order in the vicinity of the hole quasi- 
particles, i.e., the magnetic correlation length has to be 
of the order of the quasiparticle size. More precisely, our 
treatment is valid with antiferromagnetic order on length 
scales of the longest path variables included in the cal- 
culation (here 5 lattice constants). Note here that the 
magnetic correlation length in the high-T c materials is 
about 5 lattice constants for 5 = 12.5% and about 2.5 
lattice constants for 8 = 25% differing somewhat for dif- 
ferent compounds. (Numerical simulations for the t-J 
and related models yield similiar values although precise 
data cannot be obtained because of the restricted system 
sizes.) 

A remaining task would be to include doping depen- 
dent background spin fluctuations in our ansatz for the 
wave, operator Q as done in the static calculations of 
ref.E3. The scattering of quasiparticles at spin fluctua- 
tions would cause a decay of the excited quasiparticle 
states and provide lifetime broadening of the lines in the 
response functions. From this we do not expect dras- 
tic changes in the response functions for larger energies. 
A second possible improvement concerns the inclusion 
of Fermi-surface related response at low energies, work 
along this line is in progress. 



APPENDIX A: EVALUATION OF CUMULANTS 

In this appendix we show how to evaluate cumulants 
with an exponential ansatz for the wave operator Q. The 
basic relation is 



(Al) 



with At being arbitrary operators. Note that on the l.h.s. 
the operators S + , S and Ai are subject to cumulant or- 
dering whereas on the r.h.s. only the Ai operators are 
cumulant entities. However, the cumulants on the r.h.s. 
are formed with the new wavefunction = e \4>). 

Eq. ( |A1| ) can be proven either by integrating infinites- 
imal transformations (1 + S/N) and using properties of 
cumulantsEH or by explicitly using the definition of cumu- 
lant expectation values. Here we demonstrate the second 
way. We start from the definition of cumulant expecta- 
tion valuesca for a product of arbitrary operators Ai and 
an arbitrary state \<j>): 



A,=0Vi • 



(A2) 



We consider the following expression: 
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(A3) 



( = V = 
\ i =0 V i 



The last expression can be interpreted as a series expan- 
sion of the term in brackets [...] with respect to £ and r\ 
around 0: 



J3S U\c 



= n ^7 >(#^ + IIe^e 



d\i 
IP 



|Ai=0Vi 

(A4) 



In the last equation we have reintroduced (generalized) 
cumulants, now formed with the bra state {e aS (f>\ and the 
ket state \e^ s cj)). With a = (3 = 1 we obtain the desired 
result (|AL|). 



Explicitly we find from (Al) 



\e s+ Ae s \ 



\e s+ ABe s \ 



\e s+ Ae s \ 



\e s+ ABe s 



\e s+ Ae s \ 



(A5) 



\e s+ Be s 



and so on. 
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